Self-organized exotic lattices with ultracold gases 
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We study an ultracold Fermi-Fermi mixture of strongly attractive atoms trapped in a square 
lattice. We show that the interaction-induced multi-band nature of such system can change the 
inherent structure of the original lattice. The dynamically generated lattice resembles geometrically 
a Lieb lattice with interaction-driven topological insulator phases. 

PACS numbers: 67.85.Lm, 03.75.Lm, 73.43.-f 



Currently, there is a great interest in the investigation 
of topologically insulating states in solid-state physics 
P, 01 due to their potential applications in spintronics, 
topologically protected quantum computing [3], and ap- 
intomics [4]. In this regard, due to the enormous control 
and possibility of quantum engineering Q , systems of ul- 
tracold gases trapped in optical lattices provide promis- 
ing avenues to follow. One such avenue concerns optical 
creation of frustrated lattice geometries [B-Q with tun- 
able long-range hopping amplitudes, which can induce 
a variety of exotic phenomena like topological insula- 
tors, ferromagnetism etc. [8-15|. Recently, other ways 
to induce topologically insulating states such as quan- 
tum anomalous Hall states (QAH) and quantum spin 
Hall states (QSH) were investigated. Particularly in- 
teresting is the dynamical creation of interaction-driven 
topological insulators in lattices with quadratic band- 
crossmg points {QBCP) The most promising 

lattice structure to have QBCP within the Brillouin zone 
is a Lieb lattice. In [20, 21], the authors have proposed 
to create such geometry by optical means which, how- 
ever, is not easy to realize experimentally. In the spirit 
of dynamical generation of topological states, we pro- 
pose in this letter how to create such non-trivial lattice 
structures due to self-organization of the ultracold gases. 
Furthermore we show that in our system one can study 
the appearance of interaction-driven QAH and/or QSH 
states. 

We consider a mixture of two-species ultracold 
fermionic atoms trapped in an optical lattice po- 
tential Kjatt = sin^(7rx/a) -|- Va^ySin^iTry/a) + 
Vcy,zSi'n^{TTz/a), where a =t,i denotes the two species 
and Kr,a;(y)(2) the corresponding lattice depths for 
cr-fermions along the x, y, z direction respectively. The 
corresponding recoil energies are E^^ = n'^h'^ /ima-a^ , 
where denotes the mass of the a-fermion. The hll- 
ing of the tr-fermion is given by n^. For simplicity we 
consider the case in which fermionic masses are equal, 
i.e mi — m-j-. Subsequently, we define the recoil energy 
En — £'r(j. For the two-dimensional {2D) geometry we 
choose Vo = Vi^x = Vi^y, Vi = Vi^z = V^^x(y)(z), and 



Vi ^ Vq, which means that the 4,-fermions can effectively 
move in the x — y plane with the z motion freezed. The 
atoms interact with each other via s-wave scattering with 
strength as. Within the lowest-band, we can write the 
corresponding attractive Fermi-Hubbard model as Hi — 
- T.{i,j),a JoJli^ai -\U\I2 Ei,a^/ ni^ni^,. Here 
and i are the creation and annihilation operators for 
the cr-fermions in the s-band at site i = {ix,iy) with den- 
sity and the contact interaction \U\ — A'KlT?\as\/ra-c 
with mr being the reduced mass. Jqo- is the tunneling 
strength of the a-fermions in the s-band and de- 
notes nearest neighbours in the square lattice. As the 
t-fermions move in a deeper lattice (Jqj, ^ Jo^)i conse- 
quently we neglect the tunneling of the t-particles. In 
this paper we look into the spin-imbalanced situation 
with 71]^ > n-f with the excess 4,-fermions (denoted be- 
low by s-fermion) hlling m — — n-^. In the regime 
of strong attractive interaction strength, |?7| ^ Jqi, the 
'|~ and 4,-fermions will like to pair up to form compos- 



ites with creation operator 6j = sjisj; and density n^* 
22- 24 1 . Such a pair can tunnel only due to second order 
processes which are practically negligible. We can rewrite 
the Hamiltonian for the composites and s-fermions as. 
Hi = — ^oX](ij) sls^ — \U\^^nf ,wheiewehaveremoved 
the 4- subscript as only the 4,-fermions are involved in 
the dynamics. Recently it has been noted that in the 
strong interaction regime, the standard Hubbard mod- 
els have to be modified due to both intra- and inter- 



band effects |25li28| . To take these effects into account 
we write a minimal model for the excess 4,-fermions by 
including the occupation of the p-band and renormaliza- 
tion of the interaction. The inclusion of a p-band al- 
lows a 4--fermion to occupy the same site as a composite. 
The single-particle tunneling Hamiltonian then reads, 

Ht - - Jo E(ij> + Ji E(ij). piiP.j + -^i pliPyj - 

JoE{ii)yPl-Aj ~ JQY.{\i)j\\Pyr where pt,p^ (S>\.Py) 
are creation and annihilation operators of the \ particles 
in the p2:(pj,)-band. Ji is the single-particle tunneling 
amplitude in thep-band, and ^^^{y) denotes the nearest- 
neighbour sites along the x(y)-direction. Next we write 



the interaction-induced tunneling terms in the p-band as: 



E 
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(1) 



where 5,5' = x,y with 5^5' and Jn denotes intra- 
band interaction-induced tunnehng for pa;(py)-fermions 
along x{y) directions, whereas J[i denotes interaction- 
induced tunneling for p^(py)-ferniions along y{x) direc- 
tions. Such terms originate from the nearest-neighbour 
scattering due to interaction and give the most important 
contribution to the renormalization of intra-band tunnel- 
ing [29]. The interaction induced term has another com- 
ponent which can mix the s- and p-bands and can be 
written as 



Hai — -'^oi E E 

S=x,y 



.C, 



(2) 



where Jgi denotes the interaction induced inter-band 
tunneling and Cisjs — {—iy^~^^ denotes the staggered 
nature of the s — p tunneling matrix. Such hybridiza- 
tion between bands due to interaction induced tunneling 
is an important feature of the strongly interacting gases 
in optical lattices and usually has been neglected in the 
literature. Next we describe the on-site Hamiltonian for 
the excess 4,-fermions and the composites, 

Hint = -IC/2I ^ - |f/3| E ^^(^-i + ^yi) 

i i 

- \5Us\ J2 ^.i"yi^f + El + "'i)' (3) 

i i 

where n , \- = , \-V t \-- The renormalized self-energy 
of the composite is given by U2 and C/3 is the strength 
of the renormalized onsite interactions between a com- 
posite and a 4,-fermion in the p band at a site. The 
effective three-body interaction term 5U3 — U4 — 2U3 
denotes the excess energy, where C/4 is the energy of a 
single site containing one composite and one 4,-fermion 
each in the p^ and Py orbitals. In our problem we cal- 
culate the difference U3 — U2 and 5U3 within the second 
order perturbation which converges quickly with respect 
to band indices |30| and remains valid in the attractive 
interaction regime. We find that SUs is small compared 
to other parameters, so we neglect it at first. Here we like 
to comment that we derive all the tunneling and inter- 
action parameters by writing them in terms of Wannier 
functions as described in the seminal paper (3l] |. 

The total Hamiltonian H — Ht + Hqi + Hint does not 
contain any tunneling of the composites, which makes 
the commutator [hf , H] = 0. In such a case nf = 0, 1 
becomes a good quantum number for the system. The 
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FIG. 1. (Colour online) (a) The phase diagram for the con- 
figurations of the composites as a function of normalized in- 
teraction strength a = Ua/a and excess 4.-fermion filling m 
for — 1/2. The red coloured region denotes the period- 1 
checkerboard configuration for the composites. The black re- 
gion denotes the phase-separated phase, the light-blue region 
denotes the period-2 checkerboard phase, the green region de- 
notes the stripe phase and the yellow region contains mixed 
phase. The lattice parameters are Vo = 4i5H,, Vi = 30i5R.(b) 
The phase diagram with fixed m = 1/2 and composite fill- 
ing = 1/2 for Vo = 4,Er. The phases are depicted as a 
function of interaction strength a and lattice depth V\. 



ground state of the Hamiltonian H can be found by com- 
paring the energies for different configurations of nf over 
the entire lattice by solving the corresponding single par- 
ticle Hamiltonian H . We have found the lowest energy 
configurations of the composites for various parameters 
by using simulated annealing method on a 8 x 8 lattice 
with periodic boundary condition. We like to point out 
here that due to the interaction-induced origin of tun- 
neling terms in Eq.(IT]) and ©, their strengths and signs 
are defined by the various lattice depths and effective in- 
teraction strength a = as/ a. For half-filled composites 
{n^ = 1/2) the phase diagrams are shown in FigHJ In 
Figdja) we show the different phases as a function of 
m and the effective scattering length a = as /a. Vari- 
ous phases of the composites as shown in Fig. 2 are: i) 
checkerboard structure with period one ( CHI -the red re- 
gion, shown in Figl^b), ii) checkerboard structure with 
period 2 {CH2-th.e light-blue region), iii) stripe phase 
(the green region), iv) phase-separated state (the black 
region) , and v) mixed phase characterized by absence of 
any periodic structure (the yellow region). The phase- 
separated state is characterized by the clustering of the 
composites to one region of the lattice, whereas in the 
stripe phase alternative chains of the lattice are occu- 
pied by the composites. The CH2 and the stripe phase 
are only present in a small region of the phase space 
for very low filling of the excess 4,-fermions with discon- 
tinuous structural transitions. Here one should keep in 
mind that the total Hamiltonian H has similar feature 
to the celebrated Falicov-Kimball (FK) model where one 
can get a zoo of discontinuous structural phase transi- 
tions [s^l- However, normal FK models do not possess 
the multi-orbital nature and density-dependent tunnel- 
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FIG. 2. (Colour online) (a) The figure shows the distribution 
of the composites in the CHI phase. The filled (empty) cir- 
cles denotes presence(absence) of a composite, (b) The figure 
denotes the corresponding orbital degrees of freedom for the 
excess 4,-fermions. The coloured circle denotes the s-orbitals. 
The horizontal (vertical) coloured dumb-dell shapes denotes 
Pa;(py)-orbitals. (c) Pictorial description of the basis states 
for the blue (red) lattice A,B and C (A',B and C). (d) 
The ground state phase diagram corresponding to the Hamil- 
tonian ((5| as a function of dipolar strength D and contact 
interaction strength a for Vb = 4i5R,, half-filled composites 
(n^ = 1/2) and excess 4,-fermion filling m = 1/2. 

ing processes as present in our system. The mixed phase 
occurs in the region where the energy cost to occupy the 
p-orbital is small compared to other tunneling processes. 
Thus it is possible that the mixed phase is a disordered 
phase due to the composite density dependence on the 
tunneling processes. Such mixed phases have not been 



predicted before in FK like models. However one should 
be cautious about an assertion regarding the CH2, stripe 
and mixed phases in our paper, since they can be an ef- 
fect of the finite size of the system. Subsequently apart 
from the CHI and the phase-separated state, we will 
not consider other structures in the remaining part of 
this paper. In FiglTJb) we plot the phase diagram as 
a function of a and Vi with the excess 4,-fermion den- 
sity m — 1/2. The most interesting region is CHI thus 
let us mainly focus on it. The presence of CHI region 
can be qualitatively predicted when the s — p tunnel- 
ing strength Jqi becomes comparable to the inter-band 
tunneling strengths, i.e | JqiI ~ max{\Ji + Jul, Jq}, spe- 
cially for m — 1/2. The CHI phase of the composites 
is shown in Figl^Ja). The corresponding lattice struc- 
ture seen by the excess 4,-fermions is shown in the right 
panel in Figl^I^b). The coloured circular sites denotes s- 
orbitals and the dumb-bell shaped sites denote or py 
orbitals. Due to the directional nature of the inter-orbital 
tunneling Jqi in Hamiltonian ^ and the absence of any 
on-site orbital mixing term in the motion of the ex- 
cess fermions can be divided into two sub-lattices shown 
by blue and red colour in FiglSfb). Each of the red and 
blue sub-lattices have the structure of an exotic Lieb lat- 
tice. In Figm^c) we show the basis sites for the blue Lieb 
lattices denoted by A, B and C ( A',B and C for the 
red Lieb lattice ). The excess J,- fermions moving in the 
blue sub-lattice are denoted by, $i = [sA,PyB,Pxc] a-nd 
the corresponding operators in the red sub-lattice are de- 
noted by $2 — [sA'TpxBTpyc]- Their motion is governed 
by the Hamiltonian, 



H = Jf 



01 



(ij)^ {i})y {ij>x (ij>j, 

A Y (nrxi + flryi) - \SU3\ ^ 



(4) 



i,r=B,C 



i,T=B,C 



Here the first term inside the [.]-bracket in Eq.(l4]) is a 
reformulation of Hqi in Eq([2]). The second term refers to 
the energy cost of the p-orbital atoms occupying a site 
already taken by a composite, A = -Ei — IC/3I -f |J72|. The 
third term describes effective onsite interactions between 
the red and blue fermions in the sites B and C. The 
single particle dispersion relation of Lieb lattice is given 

by ek e {A, A/2 ± ^ (A/2)' + 4J2, [sin^ + sin^ ky] }, 
where momentum k = {kx,ky) belongs to reduced Bril- 
louin zone (7r/2, 7r/2). The dispersion contains Quadratic 
Band-Crossing point(QBCP) for A 5: with one of the 
dispersive bands touching the flat band at momentum 
(0,0). At A = 0, three bands touch each other at the 



momentum (0,0) with the upper and lower band having 
linear dispersion for fc <IC 1. For simplicity we consider 
only the case when A < 0. Then we can write an ef- 
fective two band Hamiltonian H = dgl + d^a^ + d^c^x, 
where (Tx(^z) a-re the Pauli matrices, X is the identity 
matrix, dp = ( Jo]^/A)(cos 2fc^ -I- cos2fcj^) and the vector 
d = {dx,dz) — {Jqi/ A){cos2kx — cos 2ky , A sin kx sin ky) . 
In this limit, the particles occupy only the B, C sites of 
the lattice and the population in the A, A' sites is negli- 
gible. For excess fermion filling m = 1/2, the dispersive 
band is filled and the system becomes a band insula- 
tor. There is a topological index W = ±2 associated 
with the vector d which makes the QBCP topologically 
stable for non-interacting systems and the corresponding 
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Berry phase is given by Wn. This shows that, due to 
interaction-induced tunnehng, one can dynamically gen- 
erate topologically protected exotic lattices from simple 
geometry. To the best of our knowledge, such proposals 
of dynamical generation of non-zero Berry phase from a 
topologically trivial non-frustrated lattice (square lattice 
in our case) have not been proposed before. 

The dynamical realization of Lieb lattice opens up an 
alternative way to study the possibility of generating in- 
teger quantum Hall effects in such systems. One possible 
way is to induce effective spin-orbit coupling, which can 
be achieved by optical means or by lattice-shaking 
[3^ . Another promising way is to generate such coupling 
dynamically by including long-range interactions which 
can create QAH and QSH states 16, Such models 
are usually hard to implement in experimentally realiz- 
able system as one needs the on-site interactions to be 
of the same order of magnitude as the long-range part 
of the interaction [l8|. To investigate such possibility in 
our system we add additional dipolar term for the excess 
4,-fermions, 



Hiuu — H + H, 



dd, 



(5) 



({iJ»,^ 



where C/dd is onsite dipolar interaction, U^y is interac- 
tion between the particles in and pj^-orbital in B and 
C sites respectively, and J7xx is next-nearest neighbour in- 
teraction between the particles in and p2:-orbital (also 
between py and py orbital) in B and C sites. 
denotes next-nearest neighbour p-orbital sites. We addi- 
tionally introduce the dimensionless dipolar interaction 
strength D — fio^'^'mi/2h?a, where /i is the magnetic 
dipole moment of the atoms and /io is the vacuum per- 
meability. In this unit, the dipole-dipole interaction has 
the form, Udd{r) — D{1 — 3z^/r^)/r^, where r is the 
inter-particle distance. For experimental realization, the 
suitable candidates are fermionic ^^^Dy which is experi- 
mentally available in quantum degenerate state [sB'l, and 
possibly fermionic ^^^Er [36]. Due to the strong attrac- 
tive contact interaction \ U\, the effect of dipolar terms on 
A is negligible. We also neglect the effective long-range 
repulsion between the composites which can further sta- 
bilize the CHI phase. Within weak-coupling limit we can 
define the mean- field parameters, 

{plBtPyC]) = -{plBzP^Cj) = iXqsH, 

(flxBi) - (flyCj) = (flxCi) - {flyBj) = XSN> 

(6) 



where Xqah denotes the order parameter for the QAH 
state that generates loop-current with broken time- 
reversal symmetry (TRS). Such state carries topologi- 
cally protected chiral-edge states. The QSH order pa- 
rameter Xqsh other hand can be thought of as two 
systems showing QAH effect with each breaking time re- 
versal symmetry, but on the whole both systems jointly 
preserve TRS. Such states contain helical edge states as 
shown in [37|. On the other hand, the spin-nematic 
state (SNS) order parameter Xsn breaks C4 symmetry 
between the blue and red sub-lattices and constitutes an 
anisotropic semi-metal [13 ■ We plot the resulting phase 
diagram in Figl2l(d) as a function of a = as/a and dipo- 
lar strength D for excess fermion filling m = 1/2. In 
the blue region the ground state is given by the spin- 
nematic, whereas in the green region the ground state 
is a topological insulator. Within the mean- field ansatz 
([6]) both QAH and QSH have the same energy, although 
this degeneracy can be broken by including higher order 
exchange interactions As the energy scale of such 
exchange term is much lower, we neglect it in our calcu- 
lations. Within mean-field theory, the temperature scale 
for the transition to the QAH / QSH state is given by, 
Tc ~ (4Jo^/A) exp[- Jo2i/2KxA] - O.OIEr for D = 0.29 
and a — —0.7. 

Let us now discuss the experimentally realizable cold 
atoms systems with the possibility of dynamically cre- 
ating exotic lattices. One such choice is the Fermionic 
^Li species in lowest two energy states with broad Fesh- 
bach resonance for the interaction strength a = as /a = 
—0.6. For the lattice constant a = 500nm results in a 
strong scattering length ~ — 300nm which is already 
achieved in Lithium mixtures in Refs. [H, Is^. For a 
mass-imbalanced mixture the effective scattering length 
is scaled and a w (as/a){l -\- m^/2m^) for the same pa- 
rameters as used in the case of equal mass. Thus, if one 
traps ''"K in the weaker lattice of Vq = AEn]^ (^-fermions) 
and ^Li in the stiffer lattice of Vi = 30£'r-|-, then the topo- 
logically protected lattice phase for a scattering length of 
as (KLi) « — 80nm can be achieved. Such a strongly at- 
tractive scattering length can be experimentally achiev- 
able in the narrow Feshbach resonance for ''''K-^Li mix- 
ture due to the tunability of the magnetic field at the 
milli-Gauss level [4o| . 

In conclusion, dynamical creation of exotic lattices can 
open up another fascinating route for experimental stud- 
ies. The proposed method is very general and can be 
extended to other lattice structures even in three dimen- 
sions. It does not involve additional optical components 
other than the ones needed for creating the parent lattice, 
which can lead to flexibility with respect to the optical 
components in experimental situation. 
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